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FUNCTIONAL CENTRAL LIMIT THEOREM FOR NEGATIVELY 
DEPENDENT HEAVY-TAILED STATIONARY INFINITELY DIVISIBLE 
PROCESSES GENERATED BY CONSERVATIVE FLOWS 

PAUL JUNG, TAKASHI OWADA, AND GENNADY SAMORODNITSKY 


Abstract. We prove a functional central limit theorem for partial sums of symmetric stationary 
long range dependent heavy tailed infinitely divisible processes with a certain type of negative 
dependence. Previously only positive dependence could be treated. The negative dependence 
involves cancellations of the Gaussian second order. This leads to new types of limiting processes 
involving stable random measures, due to heavy tails, Mittag-Leffler processes, due to long memory, 
and Brownian motions, due to the Gaussian second order cancellations. 


1. Introduction 


Let X = {X\, X 2 ,...) be a discrete time stationary stochastic process; depending on notational 
convenience we will sometimes allow the time index to extend to the entire Z. Assume that X is 
symmetric (i.e. that X = —X) and that the marginal law of X\ is in the domain of attraction of 
an a-stable law, 0 < a < 2. That is, 

(1.1) P(|Vi| > •) € RV- a at infinity; 

see iFeller (1971) or Resnick ( 1987h . Here and elsewhere in this paper we use the notation RV P for 
the set of functions of regular variation with exponent p£l. If the process satisfies a functional 
central limit theorem, then a statement of the type 


( 1 . 2 ) 


1 


|ntj 


— V x k , o < t < i 


k =1 


, 0 < t < 1 


holds, with (c n ) a positive sequence growing to infinity, and Y = \ Y(t), 0 < t < lj a non¬ 
degenerate (non-deterministic) process. The convergence is either weak convergence in the appro¬ 
priate topology on D[ 0,1] or just convergence in finite dimensional distributions. The heavy tails 
in m will necessarily affect the order of magnitude of the normalizing sequence (c n ) and the 
nature of the limiting pro cess Y. The latt er pr ocess is, under mild assumptions , self-similar, with 
stationary increments; see Lamperti ( 1962l l and Embrechts and Maeiim 3 IH2). If the process X 


is long range dependent, then both the sequence (c n ) and the limiting process Y may be affected 
by the length of the memory as well. 
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A new class of central limit theorems for longj:ange dependent stationary processes with heavy 
tails was introduced in lOwada and Samorodnitskv ( 20151 ). In that paper the process X was a 
stationary infinitely divisible process given in the form 


(1.3) 




= [ foT n (s)dM(s), 
JE 


n = 1,2, 


where M is a symmetric homogeneous infinitely divisible random measure on a measurable space 
(E,£), without a Gaussian component, with control measure /r, / : E —>• R is a measurable func¬ 
tion, and T : E —>• E a measurable map, preserving the measure /it; precise definitions of these 
and following notions are below. The regularly varying tails, in the sense of a, of the process 
X are due to the random measure M, while the long memory is due to the ergodic-theoretical 
pr oper ties of the map T, assumed to be conservative and ergodic. In the model considered in 
Owada and Samorodnitskv! ( 2015h the length of the memory could be quantified by a single pa¬ 
rameter 0 < /3 < 1 (the larger j3 is, the longer the memory). Under the crucial assumption that 


(1.4) 


M(/) == / f{s)fi(ds 
JE 


/o 


(with the integral being well defined), it turns out that the normalizing sequence (c n ) is regularly 
varying with exponent H = (3 + (1 — j3)/a, and the limiting process Y is, up to a multiplicative 
factor of the /3-Mittag-Leffler fractional symmetric a-stable (SaS) motion defined by 


(1.5) 


Y a ,p(t)=[ Mp((t- s) + ,u/)dZ a>/3 (u/,s), f > 0, 

JQ'x [0, oo ) 


where Z a p is a SaS random measure on f V x [0, oo) with control measure P 7 x vp. Here vp is 
a measure on [0, oo) given by i/p(dx) = (1 — /3)x~P dx, x > 0, and Mg is a Mittag-Lefhcr process 
defined on a probability space (IT, E l , P 7 ) (all the notions will be defined momentarily). The random 
measure Z a> p and the process Y aj p, are defined on some probability space (fl, J 7 , P). 

The /3-Mittag-Lcffler fractional SaS motion is a self-similar process with Hurst exponent H as 
above. Note that 


(1,1/a] 

He l { 1 } 

(l/a, 1) 

which is the top part of the feasible region 

(0,1/a] 
He{ (0,l] 
( 0 , 1 ) 


if 0 < a < 1, 
if a = 1, 
if 1 < a < 2, 


if 0 < a < 1, 
if a = 1, 
if 1 < a < 2 


for the Hurst exponent of a self-similar SaS process with stationary increments; see Samorodnitskv and Taqqu 
( 1994] ). This is usually associated with positive dependence both in the increments of the process 
Y itself and the original process X in the functional central limit theorem (11.211 ; the best-known 
example is that of the Fractional Brownian motion, the Gaussian self-similar process with station¬ 
ary increments. For the latter process the range of H is the interval (0,1), and positive dependence 
corresponds to the range H G (1/2,1). 

In the Gaussian case of the Fractional Brownian motion, negative dependence (0 < H < 1/2) is 
often related to “cancellations” between the observations; the statement 


J2 Cov(A 0 , X n ) = (1 


n=—oo 
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is trivially true if the process X is the increment process of the Fractional Brownian motion with 
H < 1/2, and the same is true in most of the situations in (II.2|) . when the limit process is the 
Fractional Brownian motion with H < 1/2. 

In the infinite variance case considered in Owada and Samorodnitskv ( 201 oh . “cancellations” 
appear when the integral p(f) in (ll.4jl vanishes. It is the purpose of the present paper to take 
a first step towards understanding this case, when the long memory due to the map T interacts 
with the negative dependence due to the cancellations. We use the cautious formulation above 
because with the integral p(f) vanishing, the second order behaviour of / becomes crucial, and 
in this paper we only consider a Gaussian type of second order behaviour. Furthermore, even in 
t his case our ass umptions on the space E and map T in (11.31) are more restrictive than those in 
Owada and Samorodnitskv! ( 2015l h Nonetheless, we still obtain an entirely new class of functional 


limit theorems and limiting fractional SaS motions. 

This paper is organized as follows. In Section [2] we provide the necessary background on infinitely 
divisible and stable processes and integrals, and related notions, used in this paper. In Section [3] 
we describe a new class of self-similar SaS processes with stationary increments, some of which will 
appear as limits in the functional central limit theorem proved later. Certain facts on general state 
space Markov chains, needed to define and treat the model considered in the paper, are in Section 
[H The main result of the paper is stated and proved in Section [5j Finally, Section [6] is an appendix 
containing bounds on fractional moments of infinitely divisible random variables needed elsewhere 
in the paper. 

We will use several common abbreviations throughout the paper: ss for “self-similar”, sssi for 
“self-similar, with stationary increments”, and SaS for “symmetric a-stable”. 


2. Background 


In this paper we will work with symmetric infinitely divisible processes defined as integrals of de¬ 
terministic functions with respect to homogeneous symmetric infinitely divisible random measures, 
the symmetric stable processes and measures forming a special case. Let (E,£) be a measurable 
space. Let p be a er-finite measure on E, it will be assumed to be infinite in most of the paper, 
but at the moment it is not important. Let p be a one-dimensional symmetric Levy measure, i.e. 
a (j-finite measure on M \ {0} such that 

/ min(l, x 2 ) p(dx) < oo . 

Jr 


If £q = {R £ £ : p(A) < oo}, then a homogeneous symmetric infinitely divisible random measure M 
on (E,£) with control measure p and local Levy measure p is a stochastic process (M(A), A € To) 
such that 

E e iuM(A) _ eX p | j (i _ cos ( W x)) p((ix)| mGI 

for every A € £q. The random measur e M is ind epe ndently scattered and u-additive in the usual 
sense of random measures; see Rajput and Rosihski ( 19891 ) . The random measure is symmetric 
a-stable (SaS), 0 < a < 2, if 

p(dx) = a|x| -< '“ +1 ) dx . 

If M has a control measure p and a local Levy measure p, and g : E —>• M is a measurable 
function such that 


( 2 . 1 ) 


mm 


in(l, x 2 g(s ) 2 ) p(dx) p(ds) < oo , 
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then the integral f E g dM is well defined and is a symmetric infinitely divisible random variable. In 
the a-stable case the integral is a SaS random variable and the integrability condition (12.11) reduces 
to the L a condition 

(2.2) [ |g(s)|“ g(ds) < oo . 

JE 

We remark that in the a-stable case it is common to use the a-stable version of the control measure; 
it is just a scaled version C a g of the control measure g, with C a being the a-stable tail constant 
given by 


C a = 


-l 


■-“sinx&l = (l-«)/(r(2-«)cos(W2)) 

2/vr 


if a / 1, 
if a = 1. 


See iRajput and Rosinskil (1989) for this and the subsequent properties of infinitely divisible pro¬ 
cesses and integrals. 

We will consider symmetric infinitely divisible stochastic processes (without a Gaussian compo¬ 
nent) X given in the form 

X(t)= [ g(t,s)M(ds ), teT, 

JE 

where T is a parameter space, and g(t, •) is, for each t € T, a measurable function satisfying (12.11) . 
The (function level) Levy measure of the process X is given by 

(2.3) k x = (p x g) o A'” 1 , 

with K : M x E —M r given by K(x , s) = x(g(t , s), t € T), s € A, x € R. 

An important special case for us is that of T = N and 


(2.4) g{n,s) = / oT"(s), n= 1,2,... , 

where / : E —>• R is a measurable function satisfying m, and T : E —>• E a measurable map, 
preserving the control measure g. In this case we obtain the process exhibited in (|1.3I) . It is 
elementary to check that in this case the Levy measure k x in (12.31) is invariant under the left shift 
6 on M n , 

6(x i,x 2 ,x 3 , ■■■) = (x 2 ,x 3 ,. 

In particular, the process X is, automatically, stationary. There is a close relation between certain 
ergodic-theoretical properties of the shift operator 0 with respect to the Levy measure k x (or of the 
map T with respect to the control measure g) and certain distributional properties of the stationary 
process X; we will discuss these below. 

Switching gears a bit, we now recall a crucial notion needed for the main result of this paper 
as well as for the presentation of the new class of fractional SaS noises in the next section. For 
0 < /3 < 1, let (, Sp(t )) be a /3-stable subordinator, a Levy process with increasing sample paths, 
satisfying = exp{—for 9 > 0 and t > 0. The Mittag-Leffler process is its inverse 

process given by 

(2.5) Mg(i) := Sp~(t) = infjrt > 0 : Sp(u) > t}, t > 0. 

It is a continuous process with nondecreasing sample paths. Its marginal distributions are the 
Mittag-Leffler distributions, whose Laplace transform is finite for all real values of the argument 
and is given by 

00 (9f9\ n 

Eexp «(')}=Ef^), 


( 2 . 6 ) 


6 € R; 
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see Proposition 1(a) in Bingham] ( 197ll ). Using ()2.6p . the definition of the Mittag-Lefflcr process 
can be naturally extended to the boundary cases f3 = 0 and j3 = 1. We do this by setting Mo(0) = 0 
and Mo(t) = E s t, t > 0, with E st a standard exponential random variable, and M\(t) = t, t > 0. 

The Mittag-Lefflcr process is self-similar with exponent /?. It has neither stationary nor inde¬ 
pendent increments (apart from the degenerate case /3 = 1). 

3. A NEW CLASS OF SELF-SIMILAR SaS PROCESSES WITH STATIONARY INCREMENTS 

In this section we introduce a new class of self-similar SaS processes with stationary increments. 
A subclass of these processes will appear as a weak limit in the functional central limit theorem in 
Section[5j but the entire class has intrinsic interest. Furthermore, we anticipate that other members 
of the class will appear in other limit theorems. The processes in this class are defined, up to a 
scale factor, by 3 parameters, a, /3 and 7: 

0 < a < 7 < 2 , 0 < /3 < 1 . 

We proceed with a setup similar to the one in (11.51) . Define a cr-finite measure on [0, 00 ) by 

(1 — / 3 )x _/3 dx, 0 < j3 < 1, 

5 0 (dx), /3 = 1 

(5o being the point mass at zero). Let (D',J :V ,P / ) be a probability space, and let ( S-y(t,u ')) be 
a S 7 S Levy motion and (Mp(t,ui')) be an independent ^-Mittag-Lefflcr process, both defined on 
We define 


P /3 {dx) = 


(3.1) 


Y r 


Jn' 


S~,(Mp((t — x)+,uj'),uj') dZafiiu )', x ), t > 0, 

f n'x [0,00) 

where Z a p is a SaS random measure on Q' x [0, 00 ) with control measure P / x vp\ we use the 
a-stable version of the control measure (see the remark following ( 12 . 21 ) ). 

Remark 3.1. The boundary cases /3 = 0 and (3 = 1 are somewhat special. In the case /3 = 0 we 
interpret the process in (13.11) as 

= [ S 1 (E st (uj'),oj')l(x < t)dZ a)0 (io',x), t > 0, 

J [0, cxd ) 

where E st is a standard exponential random variable defined on (fL, J 7 ', P'), independent of (5 7 (t, uj')), 
while Z a fi is a SaS random measure on Q' x [0, 00 ) with control measure P 7 x Leb. It is elementary 
to see that this process is a SaS Levy motion itself, and the dependence on 7 is only through a 
multiplicative constant, equal to 

(E / |5 7 (l)rE'(£( s “ /7 )) 1/Q . 

In the second boundary case j3 = 1 the variable x in the integral becomes redundant, and we 
interpret the process in (13.11) as 


Ya,i,-y(t) = [ S~,(t,u')dZ aA (uj'), t> 0, 
Jn' 


where Z a: 1 is a SaS random measure on with control measure PL This process is, distributionally, 
sub-stable, with an alternative representation 

Y a ,i„(t) = t> 0 , 

with W a positive strictly a/ 7 -stable random variable independenUof the SyS Lev y motion (5 7 ), 
both of which are now defined on (Cl, JyP); see Section 3.8 in ISamorodnitskv and Taqqu ( 1994] ). 
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Proposition 3.2. ( E Q „g, 7 (£ )) is a well defined H-sssi SaS process with 


(3.2) 



7 


i zA 

a 


Proof. The boundary cases /3 = 0 and ft = 1 are discussed in Remark 13.11 so we will consider now 
the case 0 < (3 < 1. 

The argument is similar to that of Theorem 3.1 in Qwada and Samorodnitskv ( 2015h . To see 
that (k" Q ,/ 3 , 7 (t)) is well defined, notice that for t > 0 , by self-similarity of (5 7 ), 


roo 

E' / |- ®)+))|“(l - fi)x-Pdx 

Jo 

roo 

= E'|S 7 (l)| a E' / Mp((t- X ) + ) ah {\ -fi)x~ p dx 

J o 

< E'|5 7 (l)|“ E< oo ; 

the finiteness of E / |S' 7 (l)| a follows since a < 7 . Next, let c > 0, ti, ...,£&> 0, and 9 \,..., 6^ € R. 
We have 



where we substituted x = cy in the last step. Because of the self-similarity of (Mg) and (S 7 ), the 
above equals 



This shows that Y a ^ a is H -ss with H given by (13.21) . 

Finally, we check stationarity of the increments of Y a ^ tl . We must check that for any s > 0, 
ti,... ,4 > 0 , and 0 i,... A € M, 



k 

Y 0j + s- x)+)) - <S , 7 (Mg((s - x)+))] 

3 = 1 


k 


Y 0 j S 'y( M l3(.( t j -*)+)) 

3 = 1 


a 

x~^dx. 


a 

x~^dx 
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Split the integral in the left-hand side according to the sign of s — x and use the substitutions 
r = s — x and — r = s — x to get 


( E' 
Jo 


+r)) — S^Mp(r))\ 

3 =1 


(s — r) 13 dr 


+ 


E' 


io 


£0AW(*i-r)+)) 

3 =1 


(s + r) 13 dr. 


Rearranging terms, we are left to check 

k 


E' 


io 


+r)) - S^Mp{r))] 


3 = 1 
k 


(s — r) 13 dr 


(3.3) 


E' 


^djS^Mpdtj -x)+)) 
j=i 


^x ^ — (s + x) ^ dx. 


However, by the stationarity of the increments of (5' 7 ), the left-hand side of (13.31) reduces to 


E' 


e 3 S ii M p{t j + r) - Mp(r)) 

3 = 1 


(s — r) 13 dr. 


Let 5 r = Sp (Mp(r)) — r be the overshoot of the level r > 0 by the /3-stable subordinator (Sp(t)) 
related to ( Mp(t )) by (12.51) . The law of 5 r is given by 

P{5 r € dx) = s ^ n r P(r + x)~ x x~^ dx, x > 0 ; 

7T 

see Exercise 5.6 in Kvprianoul ( 2006 ). Further, by the strong Markov property of the stable subor¬ 
dinator we have 

{Mp(t + r) - Mp(r ), t > 0) = {Mp((t - S r )+), t > 0) , 
with the understanding that Mp and 5 r in the right-hand side are independent. We conclude that 

k 


(3.4) 


sin /3n 


7r 


[ E' 

Jo 

roo rs 

E' 

Jo Jo 


e jS 1 {Mp(t j + r) - Mp(r )) 

3 =1 


(s — r) 13 dr 


J28 j S 1 {Mp((t j - x)+)) 

3 = 1 


r /3 (r + x) *x 13 (s — r) l3 drdx. 


Using the integration formula 


* ' t V 1 


dt = 


IT 


1 — t J t + y sin /3n 


1 - 


1 + 2 / 


, 2 / > 0 ; 


given on p. 


338 in Gradshtevn and Rvzhik ( 20041 ). shows that (|3.4j) is equivalent to (13.31) . 


□ 
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Remark 3.3. The increment process 

= YaAl {n + 1) - Y aAl {n), n = 0,1, 2,... , 

is a stationary SaS process. The argument in Theorem 3.5 in Owada and Samorodnitskv ( 2015 ) 
can be used to check that, in the case 0 < (3 < 1, this process corresponds t o a c onservative 
null operator T in the sense of (12.4jl . Furthermore, this process is mixing. See Rosiriski ( 19951 ) and 
Samorodnitskvl (120051 ) for details. On the other hand , in the case /3 = 0, the increment process is an 
i.i.d. sequence and, hence, corresponds to a dissipative operator T; recall Remark f3.il Furthermore, 
in the case (3 = 1, the increment process is sub-stable, hence corresponds to a positive operator T. 
In particular, it is not even an ergodic process. 


4. Some Markov chain theory 

The class of stationary infinitely divisible processes for which we will prov e a functional central 
limit theorem, is based on a dynamical system related to Example 5.5 in Owada and Samorodnitsky 
( 20151 ). In the present paper we allow the Markov chains involved in the construction to take values 
in a space more general than Z. 

We follow the setup of Chen ( 2001)1 ) and prove additional auxiliary results we will need in the 
sequel. Let (Z n ) be an irreducible Harris recurrent Markov chain (or simply Harris chain in the 
sequel) on state space (X, X) with transition probability P(x,A) and invariant measure n (A). Our 
general reference for such processes is iMevn aud Tweedid (2009 ). As usually , we assume that the a- 
field X is countably generated. We denote by the probability law of ( Z n ) with initial distribution 
u, and by the expectation with respect to P„. 

The collections of sets of finite, and of finite and positive vr-nreasure are denoted by 

X + := {A G X such that ir(A) > 0}, Xq~ := {A G X such that 0 < 7 r(A) < oo} . 

Since the Markov chain is Harris, for any set A G X + and any initial distribution v, on an event of 
full probability with respect to P^, the sequence of return times to A defined by ta(0) = 0 and 


(4.1) 


T^(fc) = inf{n > Tj\(k — 1) : Z n G A} for k > 1, 


is a well defined finite sequence. An alternative name for ta( 1) is simply ta- 
For a set A G X^t we denote by 


(4.2) 


n „ 

a n (v, A) = 7t(A) _1 V / P k (x, A)v(dx ), 
fc =r 7x 


n > 1 


the mean number of visits to the set starting from initial distribution u, up to time n, relative to 
its 7r-measure. When needed, we extend the domain of a to [1, 00 ) by rounding the argument down 
to the nearest integer. 

An atom a of the Markov chain is a subset of X such that P(x, ■) = P(y, •) for all x, y G a. For 
an atom the notation P n and E a makes an obvious sense. A finite union of atoms is a set of the 
form 


(4.3) 


D = U a * > q < 00 ’ 


2=1 


where cp G Xq~, j = 1,..., q are atoms. A ny suc h set is a special set, otherwise known as a D-set, 
see Definition 5.4 in Nummelin ( 1984) 1 and Orevl ( 197lh . p. 29. The importance of this fact is that 
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for any two special sets (and, hence, for any two finite unions of atoms) D± and D 2 , 
(4.4) 


a n {yi,Di) 

lim ;-—— = 1 


a n {v2,D 2 ) 


for any t wo i nitial distributions and v 2 , see Theorem 2 in Chapter 2 of Orevl (1971) or Theorem 
7.3 in fNummelinl (1984 ) (without the assumption of “speciality” there might be 7r-small exceptional 
sets of initial states). This fact allows us to use the notation at := at(v, D ) for any arbitrary fixed 
special set D and initial distribution u when only the limiting behaviour as t —>• oo of this function 
is important. For concreteness, we fix a special set D and use u(dx) = ir(D)~ 1 l£ ) (x)Tr(dx). 

A Harris chain is said to be /3-regular, 0 < /3 < 1, if the function (at) is regularly varying at 
infinity with exponent (3, i.e. 

lim act/at = c@ for any c > 0. 


£—>•00 


Let / : 


be a measurable function. For a set A E Xq the sequence 
Zk(A)= f(Zj),k = 1,2,... 


TA(fc) 

E 

j=r A (k- 1)+1 


is a well defined sequence of random variables under any law It is a sequence of i.i.d. random 
variables under V v if A is an atom and v is concentrated on A. 

The following two conditions on a function / will be imposed throughout the paper. 


(4.5) 

(4.6) 

It follows that 

(4.7) 


/ E T 1 (7t) Pi L 2 (ir) and f f(x)n(dx)=0, 

Jx 

OO 

converges in L^X, *,7r). 


k =1 


n n 

&}■■= / f 2 (x)TT(dx) + 2^ / f(x)P k f(x)ir(dx) 

Jx 4x 


< oo . 


We refer the reader to Chen ( 1999bi ) and Chenl ( 2000 ) for a discussion and examples of functions 
/ satisfying conditions (|4.5|) and (14.6h . An i mporta nt implication of the above assumptions is 
the following result, proven in Lemma 2.3 of Chen ( 1999a ): if a E is an atom such that 
inf sea \f( x )\ > 0 (i.e. an /-atom), then 


(4.8) 




We prove next a functional version of Theorem 1.3 in Chen ( 2000ll . I n infinite ergodic theory 
relate d resu lts are known as Darling-Kac theo rems; see AaronsonT i 198 il l . iThaler and Zweimiiller 
( 20061 ) . and Owada and Samorodnitskv ( 2015l l. The result below can be viewed as a mean-zero 
functional Darling-Kac theorem for Harris chains. 

Theorem 4.1. Suppose that (Z n ) is a (3-regular Harris recurrent chain with 0 < [3 < 1, and 
suppose f satisfies conditions (USD and (14.611 . Set for nt E N 

nt 

s nt (f) = Y^f(Zk), 

k =1 
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and for all other t > 0 define S n fif ) by linear interpolation. Then for any initial distribution v, 
under 


(4.9) 


(~^S nt (f),t> o) “ ((r(/3 + l)) 1/2 <J f B(Mp(t)), t > o) , 


weakly in C[0, oo), where B is a standard Brownian motion, which is independent of the Mittag- 
Leffler process Mg. If (3 = 0, then the same convergence holds in finite dimensional distributions. 


As in Chen ( 200C)ib t he proof of Theorem 14.11 proceeds in three steps: regeneration, the split 
chain method of Nummelin ( 1978h . and finally the use of a geom etri cally sampled approximation, 

Sm 

(120001 1. assuming 


also known as a resolvent approximation (see Chapter 5 in iMevn and Tweed!. 
In the first step we derive a functional version of a part of Lemma 2.3 in jCh 


npl 
ic ( 


en 


For a related result see Theorem 5.1 in Kasahara ( 19841 ). 
We define the discrete local time at a by 


existence of an atom a E Xq 


t a {ri) := max{&; > 0 : r a (k ) < n} 


Then the sequence 

(4.10) ((6(a),r a (fc) - r a (k - 1)), k = 1,2,...) is i.i.d under P a . 


Lemma 4.2. The convergence (I4.9|) holds under the assumptions of Theorem \f.l\ with the addi¬ 
tional assumption that (Z n ) has an f -atom a. 


Proof. Unless stated otherwise, all the distributional statements below are understood to be under 
P*,, for an arbitrary fixed initial distribution v. Let 

n 

H n ) = pk ( a -' a )’ n = 2 > ■ ■ ■ > 

k =1 


and note that by (14.41) . 


(4.11) 


lim <f(n)/a n = 7r(a). 
n—> oo 


Let cf 1 denote an asymptotic inverse of (j). 
stochastic process 


By Lemma 3.4 in Chen ( 1999al l. 


as n —>• oo, the 


T n (t) \= r a (\nt\)/<j) 1 (n), f>0 

converges weakly in the J \-topology on V[0, oo) to a /3-stable subordinator with the Laplace trans¬ 
form exp{— tQP/T(/3 + 1)} for 0 < /3 < 1, and to a line with slope one when /3 = 1 (we will deal with 
fi = 0 in a moment). Setting 


(4.12) 


nt 

W n {t) = n~ 1/2 ^2f, k (a), t> 0 

k =i 


(defined for fractional values of nt by linear interpolation), the laws of { {{W n (t), t > 0), ( T n (t ), t > 
0))} are tight in C[0, oo) x T>[ 0, oo) since the marginal laws converge weakly in the corresponding 
spaces as n —» oo. By (14.101) every subsequential limit is a bivariate Levy process, with one 
marginal process a Brownian motion, and the other marginal process a subordinator. By the 
Levy-Ito decomposition, the Brownian and subordinator components must be independent, so all 
subsequential limits coincide, and the entire bivariate sequence converges weakly in C[0, oo) x 
P[0, oo) to a bivariate Levy process with independent marginals. 

If 0 < j3 < 1, weak convergence of (T n (t)) is easily seen to imply, by inversion, finite dimen¬ 
sional convergence, as n —>• oo, of the sequence (£ a (nt)/<f(n)) (defined for fractional nt by linear 
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interpolation) to T(/3 + 1 )Mp. Since the pat hs are increasing, and the limit is continuous, this 
guarantees convergence in T>[ 0, oo); see Bingham : ( 197ll ) and Yamazatol (2009). Similarly, we obtain 
weak convergence to a line with slope one for 0 = 1. 

In the case 0 = 0, by Theorem 2.3 of IChenl ( 1999al ) and (|4.11j) . we see that the sequence of 
processes (£ a (nt) / 4>{n)) converges in finite-dimensional distributions to a limit, equal to zero at 
t = 0 and consisting of the same standard exponential random variable repeated for all t > 0. 
Moreover, by Lemma 2.3 in Chen (20001), the exponential random variable is independent of the 
limiting Brownian motion in (14.121) . when we perform a subsequential limit scheme for the bivariate 
process {((IFn(£), t > 0), (£ a (nt)/4>(n )))} N , similarly to the above. This time the convergence is 
in finite-dimensional distributions. 

Suppose now that 0 < 0 < 1. If (D + [0, oo) denotes the subset of D[0,oo) consisting of nonneg¬ 
ative functions, then the composition map (x, y) —x o y from C[0,oo) x D + [0,oo) to D[0,oo) is 
continuous at a point (x, y ) if y is continuous. It follows, therefore, by the continuous mapping 
theorem that 


(4.13) 


\M 


n 


la([nt\) 

■ Y £,k(a),t> 0 


k =1 


((T(/3 + l)E 0 £i(a) 2 ) 1/2 B(M^(t)), t > 0 


in D[0,oo), with (M^(t)) independent of ( B(t )) in the right hand side. By (14.111) we obtain 


1 


Y “ ((7r(a)r(/3 + l)E a ei(a) 2 ) 1 / 2 B(M /3 (t)),t>0 


k =1 


Recalling ()4.8I) . we have shown that 
(4.14) 


( , f-a{\nt\) \ 

l^= Y 6(a),t>0j n ^((T(0 + l))^a f B(M,m 


t > 0 


in D[0,oo). 

We now pr ocee d to relate (14.14j) to the statement of the lemma. First of all, Corollary 3 in 
Zweim ii llerl ( 2007l i allows us to simplify the situation and assume that the chain starts at the 
/-atom a. We can write 


nt 


l a (\nt\) 


(4.15) 


Yf( Z k)i 1 > 0 


Y Zk(a) 


k =1 


k =1 


1 — nt + [nt\ 


[nt] 

E 

k= 

TaCadnbDH 1 


f(^k) + 


nt — [nt\ 


[ntj+l 

E 

k= 

r a {(-a([nt\))+l 


\ 


f(z k ), t > 0 




Since the convergence in (14.141) occurs in the J\ topology on D[0,oo) and the limit is continuous 
(recall that we are considering the case 0 < 0 < 1), in order to prove convergence of the processes 
in (14.151) in C[0, oo), we need only show that the second and the third terms in the right hand side of 
(14.151) are negligible in C[0, oo). We treat in details the second term; the third term can be treated 
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similarly. Restricting ourselves to C[0,1], we will prove that 

[nt\ 

E /( z *> 

k= 

Ta^a([nt\)) + l 

in probability. To this end we rewrite this expression as 


sup 
/On 0<t<l 


1 


-— max 

/a n m=0,...,n 


E /(^) 


k= 

r 0 (£o(m))+l 


< 


1 


max 


' a n j=0,.../ a (n) 


max 

m= 

Ta(j)+l,-,T a (j+l) 


E /(^) 


k= 

Ta(i)+i 


Letting Wj denote the inner maximum, we must show that for any e > 0, choosing n large enough 
implies 


o i max Wj > eJUff ) < e. 

i?=0,...,4(n) 


Since the sequence (£ 0 (n)/a n ) converges weakly, it is tight, and there exists M e such that for all n, 
Pa (£a(n) > M e a n ) < e/2. Thus we need only check that for n large enough 


(4.16) 


Paf 


max Wj > e^ 
\0<j<M e a n 


= 1-P a( 


max Wj < eJaf 
\0<j<M e a n 


= 1-(1-P a (W?>e 2 a n )) 


[M e a n J + l 


< e/2. 


To see this we use Lemma 2.1 in Chen i 2000! ; which states that under our assumptions we have 
E a Wj < oo. Thus, 

Pa > e 2 On) = o(o~ 1 ) as n —>• oo , 
which verifies (|4.16l) . This proves the lemma in the case 0 < /3 < 1. 

If /? = 0, then the same argument starting with (|4.13l) works. The argument is easier in this 
case since we only need to prove convergence in finite-dimensional distributions. We omit the 
details. □ 

We are now ready to prove Theorem 14.11 in general. 


Proof of Theorem \f.l\ As before, the distributional statements below are unde rstood to be under 
P u , for an arbitrary fixed initial distribution u. The split chain method of Nummelin ( 19781 ) allows 
us to extend the result from the situation in Lemma 14.21 where we assumed the existence of an 
/-atom, to the case where we only assume that there exists a C E Xq such that 

(4.17) inf |/(x)| > 0 and P(x,A) > blc(x)TTc(A) x € X, A € X 

xeC 

for some 0 < b < 1 where vrc(-) := 7r(C') -1 7r(C' n •). 

A very brief outline of the split chain method is as follows (see Nummelinl ( 1978l l or Chapter 
5 in Mevn and Tweediel ( 20091 4 for more details). One can enlarge the probability space in order 
to obtain an extra sequence of Bernoulli random variables {Y n ). The (Y n ) are chosen so that the 
split chain (Z n , Y n ) is a Harris chain on X x {0, 1} and such that C x {1} is an /-atom (where / is 
extended to X x {0,1} in the natural way). Moreover, this can be done so that conditions (14.51) and 
(14.61) continue to hold for the split chain, and also (14.71) holds for P and n (the transition kernel 
and invariant measure for the split chain.) Then an application of Lemma 14.21 to the split chain 
proves the claim of the theorem under the assumption (14.171) . 
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The final step is to get rid of assumption (14.171) . so we no longer assume that ()4. 17|) holds to 
start with. We follow the usual procedure which, for (a small) p > 0 uses the renewal process 

N(t) := max{n > 1 : T n < t}, t > 0, 

with i.i.d. renewal intervals (r n +i — T n ), which have the geometric distribution 
(4.18) P(Id = k) := (1 - p)p k ~\ k = 1,2,... . 


The idea is to approximate the original chain (Z n ) by its resolvent chain (Zr k ), where (T^) are as 
in (14. 1 8[) . and independent of (Z n ). We then let p —>• 0. 

The resolvent chain is just (Z n ) observed at the negative binomial renewal times (T^) (this chain 
is also called a geometrically sampled chain). Its transition kernel is 

OO 

P P (x,A) := (1 -p)'^2p k ~ 1 P k {x,A), 

k =1 

and, clearly, ir is still an invariant measure. For the resolvent chain, the assumptions (14.5|) and 
(14.61) allow one to define, similarly to (14.7|) . 


a pJ 


L 


:= / f 2 (x)ir(dx) + 2 


OO 

E 

fc=i' 


f(x)Ppf(x)TT(dx) 


p OO p 

/ f 2 (x)ir{dx) + 2(1 -p) V / f(x)P k f(x)ir(dx). 
Jx ;rr [ -ix 


Furthermore, the resolvent chain is /3-regular if the original chain is, and the sequence (a[f' 1 corre¬ 
sponding to the resolvent chain (see the discussion following (|4.4|) ) satisfies 

~ (1 - p) 1_/5 a n , n ^ oo ; 


see 


(4.26) in lChen ( 200C)i ). The latter_paper also shows that the resolvent chain (Zr k ) satisfies (14. 171) 
(see also Theorem 5.2.1 in lMevn and Tweedie (2009)). 

Suppose that 0 < /3 < 1. Sinc e we h ave already proved the theorem under the assumption (14.171) . 
we can use Theorem 2.15 (c) in Jacod and Shirvaev ( 19871 ). and the “converging together lemma” 
in Proposition 3.1 of Resnick ( 20071 ) to obtain 


nt 


■E/( Z r fc ), t> 0 


k= 1 


— N(nt ), t > 0 
n 


[yj (1 - P ) 1 ” /J T(/3 + 1) a pJ B(Mp(t)), t > 0^ , ((1 -p)t, t > 0) 


in C[0, oo) x V[0 , oo). As before, it is legitimate to apply the continuous mapping theorem, to obtain 
(4.19) |--= Xj f(Z Fk ), t > 0 j (y /(1 - pY-PT(J3 + 1) a pJ B (M p (( 1 - p)t)) , t > o) 


in T>[0, oo). Repeating the same argument with the continuous version of the counting process 
(N(t)), given by 

N c (t) = N(t) + — -— for r n < t < r n+ i, n = 0,1,..., 

-L n+1 -L n 
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leads to 


N c (nt) 


(4.20) 


E f(Z Fk ), t>0j n ^° (^/ (1 - pY~P r (f3 + 1) a pJ B (Mp(( 1 - p)t)) , t> o) 


this time in C[0, oo). We now show that the convergence statement in (14.201) is sufficiently close to 
the required convergence statement in (14. 9p . and for this purpose the T>[ 0, oo) version in (14.191) will 
be useful. 

We will restrict ourselves to the interval [0,1]. Since 


— P a p,f —> <7/ as p^O, 

the second converging together theorem (see Theorem 3.5 in Resnick ( 20071 )). says that (14.91) will 
follow once we check that for any e > 0, 


(4.21) 


lim lim sup P„ I sup —— 

0 n->oo \ 0<7<1 \/Qm 


N c (nt ) 


nt 


E /( Z r,)-E/< Z *> 


> e =0. 


k= 1 k =1 

To this end we bound the probability in the left hand side of ()4.21l) by a sum of 3 probabilities: 


v i sup —L= 

o<t<i v a n 


N c (nt) 


N(nt) 


E /< z rj - E /< Z r. 


k =1 


k =1 


> 77 I +P^ I Slip — 
3 7 \ 0<7<1 V a ^ 


JV(rat) 


[ntj 


E /< z r.)-EA Z ‘) 


fc=l 


fc=l 


e 

> 3 


(4.22) 


+ Pj, | sup 
0<7<1 


fritj 


nt 


E/( z *)-E/( z *) 


fc=i 


fc=i 




Keep for a moment 0 < p < 1/2 fixed. Note that 


lim sup Pj, 

n—>oo 


sup —— 

0<t<l y/Om 


N c (nt) 


N(nt) 


E /( z r„> - E /< z r.) 


fc=l 


fc=l 


> 


< lim sup P v ( —]= max \f(Z k )\ > | ) = 0 

n—>oo \ y ^n ^^2n O y 

because, eventually, N c (n) < 2n and the convergence in (14.191) is to a continuous limit. A similar 
argument shows that for a fixed 0 < p < 1/2, 


lim sup Pj, 

71 —> OO 



iMl nt 

E/( z *)-E/( z *) 


fc=i ^=1 



= 0. 


It remains to handle the middle probability in (|4.22 jl . Let 5*, = 1 at the renewal times and 
0 otherwise. By the self-similarity of the Brownian motion and the Mittag-Leffler process, the 
convergence statement in (14.191) can be rewritten as 


NJ 


V4/(4), t >0 


™ (y/{l-p)T(fi + l)a p>f B(M fi (t)), t> 0). 













































FUNCTIONAL CENTRAL LIMIT THEOREM 


15 


Replacing p by 1 — p and, hence, each 6 k by 1 — 5 k , gives us also 

l nt i \ 

— 2(1 - 6 k )f(Z k ), t > 0 + t> 0). 

/an k= l J 

Both of these weak convergence statements take place in the J\ topology on D[0, oo). Therefore, 
N(nt) [nt] 


P sup 


0<t<l V a n 


2 f{Zr k )-^26 k f{Z k ) 


k =1 


k=\ 


> 


P sup y/pT(P + 1) (Tl-pJ \B(Mp(t))\ > - 
\0<t<l o 


which goes to 0 as p —>• 0. This completes the proof in the case 0 < ft < 1. Once again, the case 
P = 0 is similar but easier, since we are only claiming finite-dimensional weak convergence. □ 


Remark 4.3. If in Theorem 14. II the function / is supported by a finite union of atoms D, then we 
also have 

(4.23) sup E„ sup (-^=S nt (f)\ < oo 

ri> 1 0 <t<L \y/ a n J 

for the initial distribution u(dx) = ■n{D)~ l l£,{x)'K( y dx) and any 0 < L < oo. To see this, it is enough 
to consider the case where the initial distribution u is given, instead, by u(dx) = 7r(a)^ 1 l 0 (x)7r(dx), 
where a is a single atom of positive measure, forming a part of D. We use the notation in Lemma 
RTSl It is elementary to check that for each n > 1, 

£ a (n) := min{/c > 0 : r a (k) > n} = l a (n) + 1 


is a stopping time with respect to the discrete time filtration 

?k = c r (0( a )i r a(i)> j = 1,. ■ ■ ,fc), k = 1,2, ... , 


while the process 


k 

2Ci(o), A: = 1,2,... 

3 = 1 

is a martingale with respect to the same filtration. Increasing L, if necessary, to make it an integer, 
we see that 


(4.24) 


E n sup 

0 <t<L 



—E n max 

m=l,...,l a (\nL\) 



By Doob’s inequality and the optional stopping theorem, this can be further bounded by 

2 


8 

0,r 


la(\nL\) 

E n { E 0( fl )| =8E 0 (^(a)) 

3 = 1 


2 E 0 4(LnLj) 


Since 

tc ( \\2 , j E a £ a (n) 

E„(4i(a)) < oo and sup- < oo 

n> 1 tin 

by the assumption and the discussion at the beginning of the proof of Lemma 14.21 the claim (14.231) 
follows. 
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We will also need a version of Theorem EH in which the initial distribution is not fixed but, 
rather, diffuses, with n, over the set {td < n}. We will only consider the case of a finite union of 
atoms D = {JJ =l a* € <T 0 f . 

We first reformulate our Markovian setup in the language of standard infinite ergodic theory. Let 
E = X N be the path space corresponding to the Markov chain. We equip E with the usual cylindrical 
cr-field £ = A N . Let T be the left shift operator on the path space E, i.e. T(x) = (22,3:3,...) for 
x = ( 21 , 22 ,. ..) £ E. Note that T preserves the measure n on E defined by 

(4.25) /i(A) := / Pj;(-A )n(dx), for events A £ £ 

Jx 


(as usually, the notation Pj, refers to the initial distribution v = 5 X , x £ X.) Notice that the 
measure p is infinite if the invariant measure n is. This is, of course, always the case if 0 < /3 < 1. 
In the sequel we will usually assume that n is infinite even when (5 = 1. 

We will need certain ergodic-theoretical properties of the quadruple (E,£, p,T). As shown in 
Aaronson ct al. ( 1979h . T is conservative and ergodic ; this implies that 


OO 

1^4 o T n = 00 p-a..e. on E 

71=1 


for every A £ £ with p(A) > 0. For a finite union of atoms D £ Aq~ as in (14.311 . let 
(4.26) D := {x £ E : x\ £ D} 


be the set of paths which start in D. The first return time to D as defined in (14.ip can then be 
viewed as a function on the product space E via 

td ( x ) = inf{n > 0 : T n x £ D} = inf{n > 0 : x n £ D} , x = (21,22 ,...) £ E . 


The wandering rate sequence (corresponding to the set D) is the sequence h(td < n), n = 1,2,.... 
Since T is measure-preserving, this is a finite sequence. Since the Markov chain is /3-regular, this 
sequence turns out to be reg ularly vary ing as well, as we s how below. For the ergodic-theoretical 
notions used in the proof see Aaronson ( 19971 ) and Zweimtiller ( 2009 1. 


Lemma 4.4. Suppose that ( Z n ) is a /3-regidar Harris recurrent chain with 0 < /3 < 1, with an 
infinite invariant measure ir. Let D be a finite union of atoms. Then the wandering rate h{td < n) 
is a regularly varying sequence of exponent 1 — f3. More precisely, 


Ii(t d < n) ~ 


1 n 

r(i + /3)r(2-/3)^ 


as n —> 00 . 


Proof. Let Q be a Markov semigroup on X which is dual to P with respect to the measure it. That 
is, for every k = 1,2,... and every bounded measurable function / : X fc —> M, 


/ 7r(d2i) 

Jx Jx 


P(xi,dx 2 ) 


P(2fc-l,<i2fc)/(2l, . . . ,2 k) 


= / Tr(dx k ) / Q(xk, dxk-i) ■ ■ ■ / Q(x 2 ,dxi)f(xi,... ,x k ) ■ 

Jx Jx Jx 

Since ir is an invariant measure for P, it is also invariant for Q. Define a cr-finite measure on (E,£) 
analogous to p in (14.251) . but using Q instead of P, i.e. 


fi(A) := / Q x (A)7r(d2), for events A £ £ . 
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We claim that the set D given in (|4.26l) is a Darling-Kac set for the shift operator T on the space 
(E,£,p). According to the definition of a Darling-Kac set (see Chapter 3 in Aaronson ( 19971 )). it 
is enough to show that 


(4.27) 

where Tq : L\jl) 

with 


1 

“ *.l 


p(D) = tt(D) uniformly p-a.e. on D , 


L 1 (p) is the dual operator defined by 


TQ9(x) '■= 


d{flgOT 

dp, 


x 


Pg(A) := [ 5 (x)/i(dx), A € £ 

JA 

a signed measure on (E,£), absolutely continuous with respect to p. Note that the dual operator 


Tq satisfies Tgl^(x) = P k (x\,D)\ see Example 2 in i Aaronson! ([19811 ). Since we may choose a n as 
in (14.21) with A = D (recall that a finite union of atoms is a special set), it is elementary to check 
that (14.271) holds, and the uniformity of the convergence stems from the fact t hat the left side in 
(14.271) takes at most q different values on D. Applying Proposition 3.8.7 in Aaronson (1997]), we 
obtain 

I 77/ 

^ TD ~ n) ~ r(i +/?)r (2 -/3) ^' 

However, by duality, 

h(t d <n) = p(r D < n ). 

Since (a n ) is regularly varying with exponent /3, the exponent of regular variation of the wandering 
sequence is, obviously, 1 — j3. □ 


Remark 4.5. Referring to the proof of Lemma 14.41 it should be noted that using in ()4.26|) a set 
D £ Xq different from a finite union of atoms, may still define a set D that is a Darling-Kac set. 
For example, suppose that (Z^) is a random walk on R with standard Gaussian steps; that is, 


P(x, B ) = P(G £ B — x), x £ R, B Borel. 


Here G ~ A r (0,1). In this case the Lebesgue measure 7r on R is an invariant measure. It is not hard 
to see that D = [0,1] is a special set that is not a finite union of atoms. Further, a n ~ ^Jn/2^T as 
n -A- oo. We claim that the set D of paths starting in [0,1], is a Darling-Kac set for the conservative 
measure-preserving shift operator on E. To see this note that, in this case, there is no difference 
between the semigroup P and the dual semigroup Q, and, hence, 


1 n 1 n 

7^ ^ ^ klf) = 7 ^ ^ pkl[0 ’ 

v k= i v k= i 




uniformly on [0,1]. 


We can view the sums S n (f ) as being defined on the path space E by setting h(x) := f(x±), 
x = (xi, X 2 , ■ ■ ■) £ E, and then writing 

n n 

5 n (/)(x) = h ° T k ( x ) = ' 

k =1 k =1 
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This defines the notation used in Theorem 14.61 below. Define a sequence of probability measures 
On) on E by 


Vn(A) : = 


y{A n {t d < re}) 
h(t d < re) 


A e £. 


Theorem 4.6. Suppose that, in addition to the hypotheses of Theorem \f.l\ f is supported by a 
finite union of atoms D = Ui=i a * £ T’ 0 t “. Then for every L > 0, under the measures jJ. n L, 


(^=S nt {f), 0<t<Lj ™ ((T(/3 + 1 )) 1 / 2 a / 5(M /3 (t - l£)), 0 < t < L) , 

where is independent of the process B(Mp(t )) and P < x) = (x/L) 1 ^ 13 for x € [0,L] (in 
particular, = 0 a.s. if (3 = 1). The convergence is weak convergence in C[0,L] for 0 < (3 < 1 
and convergence in finite-dimensional distributions if f3 = 0. Furthermore, for all 0 < (3 < 1, 


(4.28) 


sup 

n> 1 «. 



2 

dp n < oo . 


Remark 4.7. By the similar argument to that in Remark 14.31 it is not hard to check that, under 
the assumptions of Theorem 14.61 if for some p > 2, E n |^i(D)|' p < oo for any atom a constituting 
D, then we, correspondingly, have 


sup 
n> 1 




P 

dn n < OO . 


Proof. We prove convergence of the finite-dimensional distributions first. For typographical conve¬ 
nience we will only consider one-dimensional distributions. In the case of more than one dimension, 
the argument is similar, but the notation is more cumbersome. During this proof we may and will 
modify the definition of S n t(f ) to have the sum starting at k = 0. Suppose first that L = 1 . 

Set for m = 1,2 ..., x G X and i = 1,..., q 


Pm(x,i ) := Px{Z Td = Oi\rD = m). 


Then for A £ R and a large K = 1,2,..., 
7 S nt {f) 


(4.29) 


hn 


> A 


E 


1 


—' f(td < n) 

m =0 ' 


Snt(f) 


> X,td = m) n(dx) 


\ ^ f — m) vr, 

Six ,(r D <n) £> ( *' ,)P ' 


S(nt-m)+U ) 


> A ) n(dx) 


k \kn/K\~i q 
k=l m= [(fc— l)n/K\ »=1 


S(nt—m)+ if) 


> A 


P x (t d = m) 
M ( td < n ) 


Pm(x,i ) n(dx). 


The second equality uses the fact that / is supported on D and the strong Markov property. In 
the last equality, we merely partition {0,..., n — 1} into K parts . 
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Suppose first that 0 < (3 < 1. Working backwards through an argument similar to (14.291) . we 
obtain 

K [kn/K\-l q 


E E E p ". 

k=l m=[{k—l)n/K\ i =1 


^n(t—k/K)+ (/) 


> A 


Pa: (TP = m) 
m(td < n) 


p m (x,i)Tr(dx) 


(4.30) 

where 

Clearly, 


«W/) 


— l^n 


>A , 


T ? f ,t (x) := (nt — ( nk/K — td ( x ))) + if td ( x ) € [(A: — l)n/K,kn/K). 

| nt - Tjf’^x) | 


n 


< 1/A'. 


By Theorem 14.11 and Lemma 14.41 we see that 
A' LWAI-i i 


< 4 31 > E E E p ». 

fc=l m=[(fc—l)n/Rrj i=l 


^’n(t—k/K)+ (/) 


> A 


Pa: (TD = m) 
h(td < n) 


p m (x,i)ir(dx) 


( r .. 1/9 \ /u( I (A; — l)n/A'J < td < I kn/K\ — l) 

E P (W + l)) 1/2 ^/^(A^(t - k/K)+) > A) ^- J ~ ” 1 J - L 


k=\ 

K 


E(^/Af P -({k-l)/Ky ^)p ((T(^ + ^)) l /‘ 2(T f B(Mp(t — k/K)+) > a) . 


k =1 


Combining (14.311) with (14.301) implies that 

=► (r(/8 + l)) 1 / 2 *,^^* - T 0 O ,a)+), 

Y 

where Tcxd,k is a discrete random variable independent of B and Mp such that 

P(Too,a- = k/K) = (, k/Ky-P - ((k - 1 )/Kf- p , k = l,...,K 
We claim that for every e > 0 

SU P0<5i,5 2 <t, I s l — I ^ ns l (/) Sns2 u) i 


(4.32) lim lim sup p n 

K^-og n— )-oo 


>6 = 0 . 


Then, since Too,k => as A' —>• oo, once we p ro ve (14.321) . the claim of the theorem in the case 
L = 1 will follow from Theorem 3.2 in Billingsley (119991 1. 

To see that (14.321) is true, repeat the steps in (14.291) and bound the probabilities p m (x,i) from 
above by 1. We conclude that (14.32[) is bounded from above by 

SUPo<si,S2<t, |si — S2\<.\/K | ^ns\ (/) Sns 2 (f) | X 


lim \ lim sup P ai 

K—>co ^ J n ->-oo 
1=1 


> e 


= lim q P 
K-> oo 


sup (T(/3 + VjY^crf\B(Mp(s\)) - 5(M^(s 2 ))| > 

. 0<si,S2<i, |si—S2|<l/R’ 
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where at the second step we used Theorem 14.11 Now (14.321) follows from the sample continuity of 
the process (. B(Mp(t)),t > 0). 

This proves the required convergence for L = 1. For general L we replace n by nL , t by t/L and 
use the regular variation of (a n ). Using the already considered case L = 1 we see that 


f^nL 


Snt(f) 


> A 


+ 1 )) l / 2 a f L^ 2 B(Mp(t/L - T^) + ) > A 


Since LT^ = and the process B(Mp) is /3/2-self-similar, the claim of the theorem in the case 
0 < /3 < 1 has been established. 

In the case /3 = 0 and L = 1, we proceed as in (14.291) . but stop before breaking the sum into K 
parts. Consider the case A > 0; the case A < 0 can be handled in a similar manner. Fix t > 0 and 
choose e > 0 smaller than t. Next, split the sum over m into two sums; the first over the range 
m < n(t — e), and the second over the range n{t — e) < m < nt. Denote the first sum 

g 


S„,i(A) := 




m<n(t—e) 


^’(nt—m)+ (/) 


> A 7 r(dx) 


and denote the second sum, over the range n(t — e) < m < nt. by S n , 2 (A). Let 0 < p < 1. By the 
slow variation of the sequence (o n ) there is n p such that for all n > n p and for all m < n(t — e), 
a n t-rn/ a n € (1 — p, 1 + p). By Theorem 14.11 there is h p such that for all n > n p , 


Snif) 


> 


(1±p)- 1/2 a) 


€ (1 — P, 1 + p) i 


P (a f B(E st ) > (1 ± p) _1 / 2 A) 

for each i = 1 where E st is a standard exponential random variable independent of the 

Brownian motion. For notational simplicity, we identify n p and h p . We see that for n > n p , 

u{t d < n(t - e)). 


(i -pY- 


<(1 + P) 


P (a f B(E st ) > (1 + / o)- 1/2 a) < E nil (A) 
[*{td <n) P > f 1 - ef 


h(t d < n) 

h(t d < n(t - e)) 


Furthermore, 


Letting first n 
B(E s t) that 


D n ,2(A) < 


li{n(t — e) < td < nt) 


fi(r D < n ) 

oo, then e —> 0, and, finally, p —» 0, we conclude, by the continuity of the law of 


Vn 


Snt(f) 


> A 


tP (a f B(E st ) > A) 


which is the required limit in the case /3 = 0 and L = 1. The extension to the case of a general 
L > 0 is the same as in the case 0 < /3 < 1. 

It remains to prove tightness in the case 0 < /3 < 1. We will prove tightness in C[0,1]. Since we 
are dealing with a sequence of processes starting at zero, it is enough to show that for any e > 0 
there is 5 > 0 such that for any n = 1, 2,..., 


(4.33) 


1 


p n sup 


. 0 <s,t<l, |t—s |<(5 v a » 


■]S nt (f) - S ns (f)\ >e < e. 
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However, by the tightness part of Theorem 14. 11 we can choose 8 > 0 such that for every i = 1,..., q 
and n = 1 , 2 ,..., 

_ / 1 


sup 


\0<s,t<l, |I-s|<<5 v a n 

Therefore, arguing as in (14.291) . we obtain 


\S nt (f)-S ns (f)\>e] <e. 


1 


Hn Slip 


, 0<s,t<l, |t-s|<5 V a r. 


\Snt(f)-S ns (f)\>8 


E 


P.x(td = m) 


^ ^ Pmjx, i)Pai 


1 


?n=0 


< 


V( T D < n) 

E x M(TI ,<„) Ep-(^> p - 

m=0 J A u — / j =1 


sup —— 

0<s,t<l, |t—a|<5 


| ^{nt— m)+ (/) ^(ns— m)+ (/)| >e ir(dx) 


1 


sup 


. 0<s,i<l, |i—s|<<5 V a »i 


|-?„*(/) - 5 ns (/)| > £ 7r(dx) < e. 


proving (14.331) . 

Finally, we prove (14.281) . We have 


Sn(f) 


d^n — 


1 


(4.34) 


a n li(Tn < n) J E 

1 

a n n{TD < n) 


(Sn(f)) 2 dv 
n [ f 2 (x ) 7T (dx) + 2 


n -1 n-j .. 

EE/ 

Fife," 


f{x)P k f(x) n(dx) 


where in the first step we used that / is supported by D, and in the second step we used the 
invariance of the measure n. By Lemma 14.41 and conditions (|4.5I) and (14.61) . the supremum over 
n > 1 of the right side of (|4.34[) is finite. □ 


5. A MEAN-ZERO FUNCTIONAL CLT FOR HEAVY-TAILED INFINITELY DIVISIBLE PROCESSES 

We now define precisely the class of infinitely divisible stochastic processes X = (X\. X< 2 ,...) for 
which we will prove a functional central limit theorem. Those processes are given in the form (11.31) 
of a stochastic integral. 

Let (E,£) be the path space of a Markov chain on X, as in Section |4j Let / : X — » R be a 
measurable function satisfying (|4.5|) and (14.61) . We will assume that / is supported by a finite union 
of atoms (J4.3D . Let h(x) := f(x i), x = (x\,X 2 , ■ ■ ■ ) € E be the extension of the function / to the 
path space E defined above. 

Let M be a homogeneous symmetric infinitely divisible random measure M on (E,£) with control 
measure ^ given by (14.25[) . We will assume that the local Levy measure p of M has a regularly 
varying tail with index —a, 0 < a < 2: 

(5.1) p(-, oo) € RV- a at infinity. 

Let 

(5.2) X k = [ h o T k (x) dM (x) = f /(x t )dM(x), k= 1,2,..., 

JE JE 
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where T is the left shift on the path space E. Since the function / is supported by a set of a finite 
measure, it is straightforward to check that the integrability condition ([2.111 is satisfied, so (15.211 
presents a well defined stationary symmetric infinitely divisible process. Furthermore, we have 


P(*i > A) 


/ l/( 

is 


x)\ a ir(dx) p{\, oo), A —> oo 


see e.g. Rosiriski and Samorodnitskvl ( 19931 ). That is, the heaviness of the marginal tail of the 
process X is determined by the exponent a. of regular variation in (15.1H . On the other hand, we 
will assume that the underlying Markov chain is /3-regular, 0 < (3 < 1, and we will see that the 
parameter [3 determines the length of memory in the process X. 

The main result of this work is the following theorem. Its s tatement uses the tail constant C a 
of an a-stable random variable; see ISamorodnitskv and Taqqu ( 1994 ). We also use the inverse of 
the tail of the local Levy measure defined by 


p^{y) := infjx > 0 : p(x, oo) < //}, y > 0. 


Theorem 5.1. Let 0 < a < 2 and 0 < f3 < 1. Suppose that ( Z n ) is a [3-regular Harris chain 
on (X, if) with an invariant a-finite measure n. If (3 = 1. assume that a n = o{n). Let f be a 
measurable function supported on a finite union of atoms D = u3 =1 a; E . We assume that f 
satisfies (03D and AMD - If (3 = 1, we assume also that f E L 2+£ (ir) for some e > 0. If a > 1, we 
also assume that for some e > 0, E ai \f(Z Ta2 )| 2+e < oo for any two atoms, constituting D. 

Let X = (Xi,X 2 ,...) be a stationary symmetric infinitely divisible stochastic process defined 
in g . where the local Levy measure of the symmetric homogeneous infinitely divisible random 
measure M is assumed to satisfy 40). We assume, furthermore, that 


(5.3) x P0 p(x, oo) —^ 0 as x | 0 

for some po E (0, 2). Then the sequence 

c n = C~ 1/a o 1 n /2 p^(p(T D < n) -1 ), n= 1,2,... , 

satisfies 


(5.4) c n E RVp/ 2 +(i-p)/a ■ 

Let 0 < [3 < 1. Then 

1 n 

(5.5) -^X k ^ (r(/3 + l)) 1/2 u / y aA2 (-) inC[ 0,oo), 

Cn k =l 

where (Y a ,p, 2 {t)) is the process in (13.11) . with the usual understanding that the sum in the left hand 
side is defined by linear interpolation. 

If [3 = 0, then (15.51) holds in the sense of convergence of finite-dimensional distributions. 

Proof. The fact that (15.411 holds follows from the assumption of /3-regularity and Lemma [4.41 taking 
into account that the regular variation of p at infinity implies p^~ E RV_i/ a at zero. For later use 
we also record now that 

(5.6) p(c n a“ 1/2 ,oo) r*j C a p{jD <n) 1 as n oo, 


which follows directly from the definition of the inverse and the regular variation of the tail of p in 
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We start with proving convergence of the finite-dimensional distributions. It is enough to show 
that 


1 


J ntj 


J 


=»■ (W +1)) ' 

-it-, -1 

3= 1 fc=l J=1 

for all J > 1, 0 < ti < •; • < t.i, and 8 1 ,... ,9j £ R. We use an argument similar to that in 


Owada and Samorod nitskv (|201 




The stan dard t heory of convergence in law of infinitely divisible random variables (e.g., Theorem 
15.14 in Kallenberg ( 20021 )). says that we only have to check the following: in the notation of 
Theorem 14.61 and of Section [3], for every r > 0, 


Je \ " j=i 


2 rc n\^2 9j Sntj (f)\ 1 


vp(v, oo) dv dp 


and 

(5.8) 


r 2 ~ a C a 
2 — a 


{T(/3 + l)) a/2 aJ [ [ 

J [0,oo) J £1' 


3 = 1 


P '(dui') vp(dx) 



j 

EM*,(/) 


i -1 


-l 


, oo ) dp 


r~ a C a (T(/3 + l)) a/2 cr“ 


' [0,oo) J Q,' 


3 = 1 


- z)+,u/),u/) 


P'(du>') vp(dx). 


The proof of (15.81) is very similar to that of (15.71) , so we only prove (15.71) . 

We keep r > 0 fixed for the duration of the argument. Fix also an integer L so that tj < L and 
define 


r-ry 

^{y)-=y~ 2 xp(x,oo)dx, y> 0, 

Jo 


so that the left-hand side in (15.71) can be expressed as 


A (/)!/' 

By Theorem [TU] and the Skorohod embedding theorem, there exists a probability space (fl*, J 7 *, P* 
and random variables y, y, y>, ■ ■ ■ defined on (II*, JF*,P*) such that 

-i 


p* o y- 1 = p nL o 


l 


: 'y ] OjSntj (/) J J n — 1,2,... , 


3 = 1 


p* o y- 1 = P' o | (T(J3 + l)) 1/2 <T f '£e j B{Mf ) {t j -T E ) + ) 

3 = 1 


Y n -a y, P*-a.s. 
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Then 


L 




1 OjSntj (/) | 


dp = p{t~d < nL) ( 

io* V 


\/ ®n|^n 


(IP*. 


First, we will establish convergence of the quantity inside the integral. By Karamata’s theorem 
(see e.g. Theorem 0.6 in Resnick ( 1987 )). 


(5.9) 

Therefore, as n —» oo, 
M( T D < nL) V’ 


V’(y) 




j,2—a 

2 — a 

„2—a 

2-a 

r 2—a 

2 — a 


p(y , oo) as y —>• oo. 

p(t d < nL) p(c n a“ 1,/2 |Y„|~ 1 , oo) 

|^n|"^(uD < nL) p(c n a ~ l/2 , oo), P*-a.s., 


where the last line follows from the uniform convergence of regularly varying functions of negative 
index; see e.g. Proposition 0.5 in Resnick (1987]). By (I5.6j) and the regular variation of the wandering 
rate in Lemma 14.41 we conclude that 


p(td < nL) ^ 


y/ Rn | Y n 


„2-a 


2 — a 


C a L 1_ ^|*T, P* 


-a.s. 


It is straightforward to check that 

[ L 1- ^|y|“dP* = (r(/3 + l)) a/2 cr“ 


[o,oo) J r v 


i=i 


~x) + ,w / ),a; / ) 


P ' {dkj')vp{dx), 


so it now remains to show that the convergence discussed so far can be taken under the integral 
sign. For this, we will use the Pratt lemma (see Exercise 5.4.2.4 in Resnick ( 19871 ) h The lemma 
requires us to find a sequence of measurable functions Go, G\, G 2 , ■ ■ ■ defined on (fT, J 7 *, P*) such 
that 


(5.10) 

(5.11) 

(5.12) 


p (to < nL) ^ 


< G„ P*-a.s., 


y/ n n | Y n 

G n -> G 0 P*-a.s., 
E*G n -> E*G 0 € [0, 00 ) 


Throughout the rest of the proof G is a positive constant which may change from line to line. Note 
that by (15. 6j) . p(rn < nL) i^{c n a n 7 ) tends to a positive finite constant, therefore 


p(r D < nL) i/j 


< C 




Cn^n 


- 1/2 


y n r 1 ) 


.V^\Yn\J -i/l(c n a n 1//2 ) 

Since V’ € RV- a at infinity, Potter’s bounds (see Proposition 0.8 in iResnick ( 1987] )) allow us to 


write, for 0 < £ < 2 — a\ 


ip{c n an /jl I + ) ^ > y/E^\Y n \} < C (I Y n \ a ~t + |T n r + «) 

^{c n a n ') 
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for sufficiently large n. Further, by (15.111) . y 2 ip(y) -A 0 as y | 0, which gives us 

V>( y) < C y~ 2 for all y € [0,1]. 

Thus, 


ip(c n a n 1/2 




i>(c n a n l/2 ) 
Summarizing, for sufficiently large n, 


p(td < nL ) ip 


{fin y /|| } CcL n C, 


-2 


|Tn. 


ip(c n a n 1/2 ) 


\JO'n | 


<C[ \Y n \ a -t + \Y n \ aH + 


CL n C n 


\Y n 


'Ip^CnCLn 


-!/2n 


If we we define G n to be the right-hand side of the above, then (15.10(1 is automatic. 

Let Go := (7(^1"“^ + |y|“ + £). It follows by the definition of c n and Le mm a l4~4l that 

c^a” 1 / 2 > Cp^~ ( a n /n ) —>• oo as n —> oo 

because o n /n —>■ 0 (this follows from regular variation considerations if /3 < 1, and it is assumed to 
hold if j3 = 1.) Since y 2 ip(y) —>• oo as y —>• oo by (15.9j) and the fact that a < 2, we conclude that 


CLn.C-r 


-2 


if; 


n'-n 


0 


V>(c n a n 1/2 ) 

P*-a.s., so that (15.111) holds. 

To show (15.121) . recall that by Theorem 14.61 sup n>1 E*!!^ 2 < oo. This implies uniform integra- 
bility of (|y n | a± 5, n > 1) (with respect to P*). Combining these observations, 

B*G n = C[V*\Y n \ a -t + V*\Y n \ a+ t+a n c- 2 ' 


ip(c n an 1/2 ) 


-A C(E*|y| Q -« + E*|y| a+ «) = E*<7 0 , n^oo, 

as required. This completes the proof of convergence in finite-dimensional distributions. 

It remains to prove tightness in the case 0 < f3 < 1. We start by decomposing the process X 
according to the magnitude of the Levy jumps. Denote 

pi(-) := p(- n {x : |z| > 1}) , 

P 2 {-) ■■= p{- n {x : |x| < 1}) , 

and let M t , i = 1,2 denote independent homogeneous symmetric infinitely divisible random mea¬ 
sures, with the same control measure p as M, and local Levy measures pi, i = 1,2. Then 

(W, k = 1,2,...) = (j f(x k ) ,/A/ fxj + j f(x k ) <L\J z ix), /.' = 1.2.... 

:= (xP+X™,k = 1,2,...) 

in the sense of equality of finite-dimensional distributions. Notice that and X^ 2 ) are indepen¬ 
dent. Furthermore, since / € L 2 (ir), we see that 

(5.13) E(xf } ) 2 = [ f 2 {x)n{dx) [ y 2 p(dy) < oo . 

Jx J -1 
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( 2 ) 

Fix L > 0. We will begin with proving tight ness of the normalized partial sums of X, in the 

it suffices to show that there exist 7 > 1, 


space C[0, L\. By Theorem 12.3 of Billingsley ( 1968h 
p > 0 and C > 0 such that 


(5.14) 


nt 


zx-z*, 


( 2 ) 

k 


k =1 


k =1 


A Cj) 


< £( t- s y 

- x p y ’ 


for all 0 < s < t < L, n > 1 and A > 0. 

We dispose of the case n(t — s) < 1 first, and, in the sequel, we will assume that p(td < 1) > 0. 
If this measure is zero, we will simply replace 1 by a suitable large constant 7 and dispose of the 
case n(t — s) < 7 first. It follows from (15.131) that 


nt ns 

ZX’-ZX’ 

k =1 k =1 


^ < P ^max(|x{ 2) I, |xf } |) > 


A J — 

< C\~ 2 c~ 2 n 2 (t — s) 2 


XCn 


n(t — s ) 


It follows from (15.41) that 


c n‘ r ^‘ ^ 4?V' 2 _ 2 (/3/2+(l-^)/a) — ^^1-(1-/3) (2/a-1) ■ 


Suppose first that 0 < /? < 1. If (1 — jd)(2/a. — 1) > 1, then n 2 /c 2 is bounded by a positive 
constant, and we are done. In the case of 0 < (1 — /3)(2/a — 1) < 1, since n(t — s) < 1, there is 
0 < 5 < (1 — f3){2/ot — 1) such that 

c-V (t - s) 2 < C(t - , 

which is what is needed for (15.141) . If /3 = 1, a similar argument works if one uses the stronger 
integrability assumption on / imposed in the theorem. 

Let us assume, therefore, that n(t — s) > 1. By the Levy-Ito decomposition, 


£ 4 2) 


£4 2) = / (Snt(f) - S na (f)) dM 2 

k=l JE 


= jj y{S nt (f) - Sns(f)) dN 2 + jj y{Snt(f)-S ns (f))dN 2 , 

\y(Snt (/) Sns (/))|<Ac n |?y(•S'nt (f)~Sns (/))|>Ac n 

where 1V 2 is a Poisson random measure onKxS with mean measure p 2 x p and 1V 2 := N 2 — (p 2 x p). 
Therefore, 


nt ns 

(2) \ " v (2) 
k 

k =1 /c=l 


Z^-Z 4 


< P 


+ P 


> Ac„ 


y{S n t{f) — S ns (f)) dN 2 

| y(Snt(f) — Sns (/))|^Ac n 

y{S nt (f)-S ns (f))dN 2 


^ ACr 


> 0 


|y(<Sn£ {f)—Sns (f))\>\c n 
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It follows from (15.3p that, 

P 


y(S nt (f) - Sns(f)) dN 2 

\y(Snt(f)-Sn s (f))\<Xcn 


< 


1 


A 2 4 


E 


A 2 c 2 


> A c r 


y{s nt (f) - S ns (f )) dN 2 

\y(Snt ( f)~Sns (/))|^^ c n 

[y{S nt (f)-S ns (f))] 2 dp 2 dp 


< 4 


| y{S n t(f) — Sns (/)) | ^ ACtt, 


Snt(f) - Sns(f) 
Ac,, 


/ XCn/\S nt (f)-S ns (f)\ 

/ yP 2 {y,oo)dy I dp 

V 


< 




Similarly, 


y{s nt (f) - 5 ns (/)) div 2 

| y (Snt {f) — Sns (/))|> Acn 


> 0 


< p 


’«*(/)-5 na (/))| >Xc n }) >1 


< EIV 2 ({|y(5 nt (/) - S ns (f))\ > Ac n }) 


= 2 


< 


[ P 2 {^c n \S nt (f) - S ns (f))\ 1 , oo) dp 
Je 

« [ \Snt(f)-S ns (f))rdfi. 

' Je 


C 1 


A Po dn JE 

Elementary manipulations of the linear interpolation of the sums and (14.281) show that 

C 1 


nt ns 

E4 2> -E4 2> 


> A c n < - 5 T- max 

I \P° Cn n(t—s)—l<m<n(t— s)+l, 


[ \S m (f)\ Po dp 

k=i k=i / • - Je 

< (, a n(t-s )) P ° /2 P( T D < n(t - S)) , 

where at the last step we have used the assumption n(t — s) > 1, regular variation, and Theorem 


Suppose first that 0 < (3 < 1. Choose e > 0 so that 2/a — e — 1 > 0. Note that, if (15.31) holds for 
some po E (0,2), it also holds for all larger pq. Thus, we may assume that po is close enough to 2 
so that 


( 1 -/ 3 ) 


V o ep o 


a 2 

Since we are assuming that p(td < 1) > 0, we see that 


-1 >0 1 - 


P o 


1 


-jo ( a n{t _ s) y o/2 p{r D < n(t - s)) 


0 

{ -"n. 
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< c 


< c 


f a n p(T D < n) 2/c 

\ 4 


Po /2 


( t - s )\ po/2 ( p{r D < n{t-s )) 


p 0 /a-ep 0 /2 


Kt d < n) 


(i t-s )\ P0/2 (K t d < n{t-s )) 


po/a-epo/2 


p{j D < n ) 


The first inequality above uses the choice of e and j)q. while the second inequality follows from the 
definition of c n and regular variation of . Next, we choose 0 < £ < min{/3,1-/3} such that 

'Po m 


( 1 -/ 3 ) 


a 2 

By the regular variation of a n and p(td < n), 


?)>«■ 


A*(rg<n(* a)) On^-si < c{ t_ s) p-t 

p(td <n) 


(In 


Combining these inequalities together, we have 


nt 


£ 4 2) -£a 


( 2 ) 


> Ac. 


n , < 


c_ 

A po 


(t - 8)\ 


k =1 fc=l 

where 7 = (/3 — £)po/2 + (1 — (5 — £)(po/ct — epo/2). Due to the constraints in e, po, and £, it is 

( 2 ) 

easy to check that 7 > 1. This establishes tightness for the normalized partial sums of X, in the 
case 0 < (3 < 1. 

If 0 = l,_then the assumption a n = o(n) and a standard modification of Theorem 12.3 of 
Billingsley ( 19681 ) make the same argument go through. 

It remains to prove tightness of the normalized partial sums of in the space C[0,L], for a 
fixed L > 0. For notational simplicity we take L = 1. We will consider first the case 0 < a < 1. 
Let p^~{y) := inf '{x > 0 : pi(x, 00 ) < ?/}■ y > 0 be.the inverse of the tail of p\. We will make use of 
a certain series representation; see Rosiriski ( 1990 1: 


(5.15) 


itx\ o < * < 1) = (fl wT 0 <«< 1) . 


where (< 7 ) is an i.i.d. sequence of Rademacher variables (taking +1, —1 with probability 1/2), Tj is 
the jth jump time of a unit rate Poisson process, and (vj n ' > ) is a sequence of i.i.d. random variables 
with common law p n . Further, (< 7 ), (F^), and (vj n ^) are taken to be independent with each other. 
Fix £ E (0,1 /a — 1) and for K > 2(1 /a + £) — 1, we split the right-hand side above as follows. 

K 


r( 


3 = 1 


\MX1X )) S ”‘ (/)( h 




T$(t)= J2 e iP" 




j=K +1 

We will prove that the sequence (cjpTpp) is, for every K , tight in C[0,1], while 






(5.16) 


lim limsupPf sup \T^ 2 (t)\ > ec n ) = 0 , for every e > 0 . 

K—>00 n—>00 0 <t<l 
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Notice that 




3 = 1 


\2h{t d < n) 


rww<»r , ) ^^ 


Since and are both regularly varying at zero with exponent — 1/a, 


Pi 




n —>• oo , a.s.. 


On the other hand, by Theorem 14.61 each a n 1 ^ 2 S n t{f){V^) weakly converges in C[0,1], and thus, 

by independence, c ~ 1 ±4 turns out to be tight in C[ 0,1]. 

Next, we will turn to proving (15.161) . The probability in (15.161) can be estimated from above by 

P j>r /"-W’B S ")- 1 ) o S s Up|a- 1/ ^ n «(/)(^ , )| > 

Appealing to Potter’s bounds and the fact that p\ has no mass on {x : |x| < 1}, 

" r ( 2 p(r l‘< „> ) s »r‘) < ■ 


Combining this bound, Chebyshev’s inequality, and the Cauchy-Schwarz inequality, the probability 
in (15.161) is bounded from above by 

2 


Ce~ 2 E(B n ) 2 E Y max {r 




yj = K+l 

where B n = sup 0<t<1 \a n 1 ^' 2 S n t(f)(vY)\ ■ We know from Remark 14.31 that the sequence (E(R n ) 2 
is uniformly bounded in n. Because of the restriction in K, 

2 / \ 2 
2(l/a±Ol 1/2 


E £ r 


-(l/a±0 

3 


< 


T 


d=K+l 


<c[ y r (1/Q±0 

\3 = K+\ 


£ ( Er 7 

j=K +1 

where the rightmost term vanishes as K —>• oo, so the proof of the tightness has been completed. 

This proves tightness in the case 0 < a < 1, and we proceed now to show tightness of the 
normalized partial sums of At 1 in the space C[0,1], for the case 1 < a < 2. Recall that, in this 
case, we impose a stronger integrability assumption on /. We start with the Levy-Ito decomposition 
(1 5.15 [) and write, for K > 1, 


nt 


nt 


Y 4 1} , o < t < i u ( y + 4 2 ’ K) ) ,o<t<i 

\k=1 ) 


\k=\ 


\ 


fySntif) dN\ + JJySnt(f) dNi, 0 < t < 1 

I dn | V | >KCn /CLn J 

Note that the probability that the process (X]fcf=i 4? K 0 does not identically vanish on the interval 
[0,1] does not exceed 

P(W{(x, y) : \y\ > Kc n /y/a t d (x) < n} > l) 
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< EJVi{(x, y) : \y\ > Kc n /y/a. r D (x) < n} 
= 2 / 0 i (Kcn/y/a^, oo) n(r D < n) 

P (. Kc n /yfan, oo ) 


2 C a 


2C a K~ 


P (cn/y /®nj OO ) 

as n —>• oo, and this can be made arbitrarily small by choosing K large. Therefore, we only need 
to show tightness, for every fixed K , of the normalized partial sums of the process X^' K \ As in 
(15.141) . it is enough to prove that there exist 7 > 1, p > 0 and C > 0 such that 


(5.17) 


nt ns 

E4 hK> 

k =1 k =1 


> Ac n j < -^(f-s) 7 


for allO<s<f<l,n>l and A > 0. In a manner, similar to the one we employed while proving 
(15.141) . we can dispose of the case n(t — s) < 1, so we will look at the case n(t — s) > 1. 

Let 0 < e < 1 be such that / € L 2+£ (-7r). By Proposition 16.21 

«« 2-{-e 


(5.18) 


1 


< 


+ 


cl +£ 

C 

C 2 n +€ 

c 

C 2 n +€ 


E 


nt 




(1,K) 

k 


k =1 


fc=l 


f | S n t (/) - <5 ns (/)| 2+ "|y| 2+£ l { |j / |<^ Cn/ y^ } d / 9id^ 

JRxE 

[ \Snt(f ) - 5'ns(/)| 2 |y| 2 l{| ? /|<A'c rl /ya7}rfPl^7 t 
JlxE 


l+e/2 


By Karamata’s theorem, 


\y\ 2+£ h\y\<Kc n /V^} d Pi ^ C ^T D 7 n ) ■ 


2+6 


'n / V 


Further, by the fact that / is supported on ZD, and by the integrability assumption on /, we know 
that 

/ | <Snt (/) - <S r ns(/)| 2+£ d/z < C/u(t£) < n(t - s))(a n(t _ s) ) (2+e)/2 ; 

J E 

see Remark 14.71 By Lemma 14.41 p(r£> < n)an + is regularly varying with exponent bigger than 
1, so the first term in the right hand side of (|5.18p is bounded from above by C(t — s) 7 , for some 
7 > 1. A similar argument produces the same bound for second term in the right hand side of 
(15.181) . Now an appeal to Markov’s inequality proves (15.171) . □ 

6. Appendix: Fractional moments of infinitely divisible random variable 

In this appendix we present explicit bounds on the fractional moments of certain infinitely divis¬ 
ible random variables in terms of moments of their Levy measures. These estimates are needed for 
the proof of Theorem 15.II We have not been able to find such bounds in the literature. Combinato- 
rial identities for the integer moments have been known at least since B assanand and Bona ( 199C L 
Fractional moments have been investigated, using fractional calculus, bv IMatsui and Pawla s (2014), 
but the latter paper does not give general explicit bounds of the type we need. 

We will consider fractional moments of nonnegative infinitely divisible random variables and of 
symmetric infinitely divisible random variables in the ranges needed in the present paper, but our 
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approach can be extended to moments of all orders. We start with nonnegative infinitely divisible 
random variables with Laplace transform of the form 


( 6 . 1 ) 


Ee~ ex = exp <* - 


|~ (l - e u+(dy)| := e I[e \ 9 > 0, 


with the Levy measure satisfying 


f 


yv+{dy) < oo. 


Proposition 6.1. Let 1 < p < 2. Then there is c p € (0, oo), depending only on p, such that for 
any infinitely divisible random variable X satisfying (16.11) . 

p\ 


( 6 . 2 ) 


EX P < c p 


y p v+(dy) + 


V u +(dy) 


Proof. If the pth moment of the Levy measure, 

V p v+{dy) , 

is infinite, then so is the left hand side of m, and the latter trivially holds. Therefore, we will 
assume for the duration of the proof that the pth moment of the Levy measure is finite. We reserve 
the notation c p for a generic finite positive constant (that may depend only on p), and that may 
change from line to line. We start with an elementary observation: there is c p such that for any 
x > 0, 

POO 

x p = c p (l — e~ xy ) 2 t/ -( p +1) dy . 

Jo 

Therefore, 

P OO P OO 

(6.3) EX p = c p J E (1 - e~ yX ) 2 y~ {p+l) dy = c p J (l - 2e~ I{y) + e" /(2y) ) y~ {p+1) dy , 

where I is defined in (16.11) . Denote 

9 + = sup{0 > 0 : 1(9) < l} € (0, oo] . 



Observe that 

(6.4) 9 + > Ijv+(dy) 

To see that, notice that, if 9 + < oo, then 

POO 

1 = I(9 + )<9 + yv + (dy). 
Jo 

We now split the integral in (16.31) and write 



EX p = 




:= A + B. 


B < c p 


16+ 


y 


-(p+i) 


dy < c p 


yv+(dy) 


Note that by (16.41) . 
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Next, using the inequality 1 — e 29 < 2(1 — e 9 ) for any 9 > 0, see that 1(9) < 1(20) < 2 1(0) for 
each 9 > 0. Note also that for 0 < b < 2a we have 


1 - 2e -a + e~ b < 2 a 2 + (2a - b), 


and we conclude that 


r e+ 


\ 2 pQ+ / noo \ 

(l - e~ xy ) u + (dx)j y _(p+1) dy+c p J (j (l - e~ xy ) 2 v+(dx)J y~ (p+1) dy . 


A < c p 

Jo \J 0 

Using the fact that for 0 < y < 9 + we have 


I(y) < min ^l,y f xv + (dx)^j , 


it follows that 
r e+ 


/ o 


(l — e xy ) z '+(dx)^J y {p+l) dy < J 

+ [ 

J( 


(fo°° xv+(dx)) 


y 2 [ / xv+(dx)] y (p+1) dy 


Finally, 


(So° x "+( dx )) 

r 

/ 0 


y ^ p+1) dy<c p 


poo 

J 0 


< 


j™ (1 - e~ xy f u + (dx)^j y- (p+1) dy 

poo / pOO \ /*CXD 

l (l~e~ xy ) y~ [p+1) dyj v + (dx) = c p y p v+(dy ), 


and the proof is complete. 


□ 


We consider next a symmetric infinitely divisible random variable, with characteristic function 
of the form 



for some symmetric Levy measure u, satisfying 



v(dy) < oo . 


Proposition 6.2. Let 2 < p < 4. Then there is c p € (0, oo), depending only on p, such that for 
any symmetric infinitely divisible random variable Y satisfying (|6.5D . 


( 6 . 6 ) 


E|17 < c P 


( r \y\ p v(dy) + 
\J—oo 



Proof. Once again, we may and will assume that the moments of the Levy measure in the right 
hand side of (|6.6D are finite. We start with the case when u(R) < oo. If (Wj) is a sequence of i.i.d. 
random variables with the common law u/u(M), independent of a Poisson random variable N with 
mean i/(R), then 

N 

3 = 1 
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and so by the Marcinkiewicz-Zygmund inequality (see e.g. (2.2), p. 227 in Gut (.2009)), 

p/2 


N 


E|Y| P < CpE Wj 


d = 1 


The random variable 


N 


* = 5>i 


3 = 1 


is a nonnegative random variable with Laplace transform of the form (16.11) . with Levy measure v + 
given by 

v+(A) = v{y : y 2 £ A}, A Borel. 

Applying Proposition [64] (with p/2), proves (16.61) in the compound Poisson case z/(M) < oo. In the 
general case we use an approximation procedure. For m = 1,2,... let u m be the restriction of the 
Levy measure v to the set {y : \y\ > 1/m}. Then each v m is a finite symmetric measure. If Y m 
is an infinitely divisible random variable with the characteristic function given by (16.51) . with v m 
replacing u, then Y m 4 7 as m -> oo, and the fact that (16.61) holds for Y follows from the fact 
that it holds for each Y m and Fatou’s lemma. □ 
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